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Abstract 

Current astronomical measurements indicate that approximately 73% of the uni- 
verse is made up of dark energy. Stochastically quantized self-interacting scalar 
fields can serve as suitable models to generate dark energy. We study a particular 
model where the scalar field theory underlying dark energy exhibits strongest pos- 
sible chaotic behaviour. It is shown that the fluctuating chaotic field momenta obey 
a generalized (nonextensive) statistical mechanics, with an entropic index q given 
by q = 3, respectively q = — 1 if the escort formalism is used. 



1 Introduction 

The formalism of nonextensive statistical mechanics has been developed as 
a useful tool to statistically describe complex systems [1,2,3]. The concept 
is based on generalizing the Shannon entropy Si to more general informa- 
tion measures, such as the Tsallis entropies S q [l], or other generalized in- 
formation measures as well[4,5,6], and doing statistical mechanics based on 
these more general information measures. Often, the formalism is relevant to 
nonequilibrium systems with spatio-temporal parameter fluctuations, which 
exhibit a superposition of different statistics on different time scales (in short, 
a 'superstatistics') [6,7,8,9]. In recent years, several physical applications have 
been pointed out for the formalism with q ^ 1, among them being turbu- 
lent flows[8,9], pattern forming systems[10], heavy ion collisions [11], cosmic 
rays[12,13,4], cosmology [14], econophysics[15], and biological systems [16]. 

In this paper we want to point out yet another potential application of the 
nonextensive formalism, namely to stochastically quantized dark energy mod- 
els. There is growing observational evidence that the universe is presently 
dominated by some unknown form of vacuum energy with equation of state 
close to w — — 1, so-called dark energy. Current astronomical data indicate 
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that there is approximately 73 % dark energy, 23 % dark matter, and 4% 
ordinary matter in our universe [17, 18]. The nature and origin of the dominat- 
ing dark energy component is not understood, and many different models are 
currently being discussed [19,20,21]. 

The simplest way to generate dark energy is via a self-interacting scalar field. 
The dark energy density is then essentially given by the potential of the field, 
in either a static way as for a true cosmological constant or in a dynamical 
way as in quintessence models[19,20]. In this paper we consider a model for 
dark energy which associates dark energy with self-interacting scalar fields 
corresponding to a <£> 4 -theory, which is second quantized [21,22,23]. The differ- 
ence to other scalar field approaches is that our fields are very strongly (rather 
than weakly) self-interacting, and that they are second-quantized. We use as 
the relevant method to quantize the scalar fields the stochastic quantization 
method introduced by Parisi and Wu [24,25]. In the fictitious time variable of 
this approach, the fields behave in a deterministic chaotic way. Our physical 
interpretation is to associate the chaotic behaviour of the scalar fields with 
vacuum fluctuations [22,23]. When quantum mechanical averages are formed, 
the fields generate quite smoothly distributed dark energy [21]. 

We will show that the probability density of the fluctuating chaotic momenta 
in our model is given by the generalized canonical distributions in the for- 
malism of nonextensive statistical mechanics. The relevant entropic index is 
q = 3 (see also [14]), or q = — 1 if the escort formalism is used. This means 
the dark energy component behaves quite different from a system of particles 
described by ordinary statistical mechanics. The expectation of the scalar field 
potential plays the role of a generalized thermodynamic potential in this set- 
ting. We will show that a theory formulated in terms of escort distributions, 
i.e. q = —1, has many advantages as compared to q — 3. We will explicitly 
calculate quantities like the entropy S q , the generalized internal energy U q and 
the generalized free energy F q for our model. 

If our model is indeed the correct description of dark energy in the universe, 
then this means that Boltzmann-Gibbs type of statistical mechanics (q — 1) 
is only relevant for a minority of the contents of the universe (ordinary and 
dark matter, together 27 %), whereas the dominating dark energy contents of 
the universe (73 %) is described by a different type of statistical mechanics, 
with q — — 1. 

This paper is organized as follows. In section 2 we first recall how to second- 
quantized a scalar field dynamics using stochastic quantization. We then show 
that a chaotic dynamics can be obtained if the potential is very strong. Our 
main example in section 3 is a <^ 4 -theory leading to an effective dynamics given 
by weakly coupled 3rd order Tchebyscheff maps. This dark energy model ex- 
hibits strongest possible chaotic behavior. In section 4 we show how the pa- 
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rameters have to be chosen in order to reproduce the currently measured dark 
energy. Finally, in section 5 we embed the chaotic field model into the formal- 
ism of nonextensive statistical mechanics and calculate various thermodynamic 
quantities. 



2 Stochastic quantization of strongly self-interacting scalar fields 



Let us consider a homogeneous self-interacting scalar field f in Robertson- 
Walker metric. To second- quantize it, one can use the Parisi-Wu approach 
of stochastic quantization [24,25]. This means one considers the following 
stochastic differential equation: 

^- s <p = (p + ZH<p + V'{ip) + L(t, s) (1) 



Here H is the Hubble parameter, V is the potential under consideration, t is 
physical time, and s is an artificial coordinate (called fictitious time) which 
is just introduced for quantization purposes. L(t, s) is Gaussian white noise, 
5-correlated both in t and s. The expectations with respect to the stochas- 
tic process generated by this stochastic differential equation for s — > oo are 
identical to the quantum mechanical expectations of the field theory under 
consideration. 

We may discretize eq. (1) using s = nr and t = i8, where n and % are integers 
and r is a fictitious time lattice constant, 5 is a physical time lattice constant. 
The continuum limit requires r — * 0, 5 — > 0, but we will later see that it 
makes physical sense to keep small but finite lattice constants of the order of 
the Planck length. We obtain 

^ = - 2 ^ + + 3 f " V^ 1 ) + + noise{2) 



equivalent to 



+ t ■ noise, (3) 

where a dimensionless coupling constant a is introduced as a := |j. We also 
introduce a dimensionless field variable by writing ip l n = $ l n p max , where 
p m ax is some (so far) arbitrary energy scale. The discretized stochastically 
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quantized system is equivalent to a coupled map lattice of the form 

K + i = (1 " a)T(0 + \H5a(V n - *™) + + S^ 1 ) + r ■ noise, (4) 



where the local map T is given by 

T($) = $ + ^ rF'(p mal f). (5) 



Here the prime means ' — ■§- — ]T~~M>- ^ or l ar S e ^ universes) the term 
proportional to H can be neglected and one obtains a symmetric diffusively 
coupled map lattice of the form 

¥ n+1 = (1 - a)T(¥ n ) + + S^ 1 ) + r no* S e. (6) 



Apparently the iteration of a coupled map lattice of the form (6) with a given 
map T has physical meaning: It means that one is considering the second- 
quantized dynamics of a self-interacting real scalar field if with a force V 
given by 

V'(<p) = — \-<p + Pmax T (^-) ) , (7) 

T { \Pmax J I 



equaivalent to 

V(ip) = J -\p 2 +p m ax J dipT i—^—J | + const. (8) 



3 Chaotic (//-theory 



It is interesting to see that our model can exhibit chaotic behaviour for very 
strong forces V ~ r _1 . As a distinguished example, let us consider the map 

$ n+1 = T_ 3 ($„) = -4$3+3$ n (9) 



on the interval $ G [—1,1]. T_ 3 is the negative third-order Tchebyscheff map, 
a strongly chaotic map conjugated to a Bernoulli shift of 3 symbols. This map 
arises in the scalar field dynamics (6) if the scalar field potential is given by 

V- 3 (<p) = — L 2 - -^-A + const, (10) 

^~ t Pmax J 
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or, in terms of the dimensionless field $, 

V-M = —pL^ 2 ~ $ 4 ) + const. (11) 



Apparently, this potential describes a field ip that rapidly fluctuates in ficti- 
tious time on some finite interval . The small noise term in eq. (6) can be 
neglected as compared to the deterministic chaotic fluctuations of the field. 

Of physical relevance are the expectations of suitable observables with respect 
to the ergodic chaotic dynamics. For example, one can calculate the expec- 
tation of the potential V_ 3 , which describes the vacuum energy generated by 
this chaotic field theory: 

1 — n 

<V-s{<p)) = — PL,«* 2 > - <$ 4 » + const (12) 



For uncoupled Tchebyscheff maps (a — 0), expectations of any observable A 
can be evaluated as the ergodic average 

+i 

(A) = [ (13) 



with the natural invariant density being given by 



**> = whw- (14) 



From eq. (14) one obtains ($ 2 ) = \ and ($ 4 ) = |, thus 

(V-sCp)) = \— + const. (15) 

O T 



Alternatively, we may consider the positive Tchebyscheff map T 3 ($) = 4$ 3 — 
3$. This basically exhibits the same dynamics as T_ 3 , up to a sign. Repeating 
the same calculation we obtain 

V^) = — pL,(-2$ 2 + $ 4 ). (16) 



For the expectation of the vacuum energy one gets 
1 — n 

<Vs(<f)) = —pL*(-2<* 2 ) + (^)) + const, (17) 
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which for a = reduces to 

(VsM) = 4— + const. (18) 



Symmetry considerations between T_ 3 and T 3 suggest to take the additive 
constant const as const = +^-p 2 ma x\{^ 2 ) ■ m this case one obtains the fully 
symmetric equation 

{v±M) = ^~P 2 max {-\(&) + <$ 4 )} , (19) 



which for a — > reduces to 

(^)) = ±% (4) ■ (20) 



The above vacuum energy changes in a nontrivial way with the coupling a. 
Extensive numerical simulations were performed in [22,23], and it was found 
that as a function of a it has several local minima that numerically coincide 
with running standard model coupling constants, evaluated at the mass scales 
of known fermions and bosons. This emphasizes the physical relevance of the 
model. The role of chaotic fields in the universe could be to fix and stabilize 
fundamental constants as local minima in the dark energy landscape. 



4 Constructing a dark energy model 



Let us now fix the free parameters in our chaotic field model in such a way 
that the currently experimentally measured properties of dark energy come 
out correctly. 

Firstly, the measured dark energy density is positive, not negative. This means 
we have to choose the negative Tchebyscheff map T_ 3 . 

Next, consider the parameter r. It is the lattice constant of fictitious time s and 
has dimension GeV~ 2 . Ordinary stochastic quantization based on Gaussian 
white noise requires the continuum limit r — > 0. But since quantum field 
theory runs into difficulties at the Planck scale m P i and is expected to be 
replaced by a more advanced theory at this scale, it is most reasonable to take 
the small but finite value 

t ~ m P ] ~ 10- 38 GeV- 2 . (21) 
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Next, the classical equation of state of dark energy should be close to w — — 1, 
as experimentally measured [17,18]. The classical equation of state is defined 
by w = (p) I (p), where p is the pressure and p the energy density. As shown in 
[21], this condition requires that the Tchebyscheff maps are weakly coupled. 
a = generates an equation of state exactly given by mj = —1, and a « 1 
generates an equation of state w ~ — 1. 

Finally, the parameter p max can be determined from the currently measured 
dark energy density [17,18] 

p° bs = (2.9 ± 0.3) • 10- 47 GeV 4 . (22) 

For a«0, eq. (20), (21) and (22) fix the parameter p max as 

Pmax ~ W' 42 GeV. (23) 



In static models, where the vacuum energy density does not change with the 
expansion of the universe, the smallness of this parameter is the famous cos- 
mological constant problem [19,20]. In models with dynamically evolving dark 
energy, the (current) smallness of this parameter is related to the fact that 
the universe is old [21]. 



5 Generalized statistical mechanics 



The type of chaotic vacuum fluctuations considered here have rather interest- 
ing properties, which can be described in the language of statistical mechanics. 
For a = the natural invariant measure describing the probability distribu- 
tion of the fluctuating variables $^ is given by eq. (14), and the measure 
factorizes at all lattice sites. The relevant probability density can be regarded 
as a generalized canonical distribution in non-extensive statistical mechanics 
[1]. As it is well known, one defines for a dimensionless continuous random 
variable X with probability density p(x) the Tsallis entropies as 

S q = -L [ (i-J p ( x )*dx). (24) 



Here q is the entropic index. The Tsallis entropies contain the Shannon entropy 
Si as a special case for q—>l. Extremizing S g subject to the constraint 

J p(x)E(x)dx = U (25) 
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one ends up with g-generalized canonical distributions. These are given by 
p(x) ~ (l + (q- l)(3E(x)y^, (26) 



where E is the energy associated with microstate x, and (3 is the inverse 
temperature. Of course, for q — > 1 one obtains the usual Boltzmann factor 



e-P E . 



Alternatively, one can work with the escort distributions, defined by [26] 

™ - Mr, (27) 



If the energy constraint (25) is implemented using the escort distributions P(x) 
rather than the original distribution p(x), one obtains generalized canonical 
distributions of the form 

P(x)~(l + (q-l)PE(x))-£i. (28) 



Again the limit q — > 1 yields ordinary Boltzmann factors e l3E . 

Let us now apply this formalism to the chaotic fields X = We might 
identify E = |m$ 2 as a formal kinetic energy associated with the chaotic 
fields. We then get by comparing eq. (14) and (26) 



q = 3 (29) 
p- 1 = -m (30) 

Two problems arise: Firstly, the temperature that formally comes out of this 
approach is negative. Secondly, and more seriously, the Tsallis entropy of the 
distribution (14) as defined by the integral eq. (24) does not exist, since the 
integral /^(l — x 2 )~ 3 / 2 dx diverges. 

However, it is remarkable that these problems do not occur if the escort for- 
malism [2] is used. By comparing eq. (14) and (28) we obtain 



q = -l (31) 
(3~ l =m. (32) 

We obtain the result that our dark energy model behaves similar to an ideal 
gas in the nonextensive formalism but with entropic index q — — 1, rather 
than q = 1 as in ordinary statistical mechanics. The 'velocity' v is given by 
the dimensionless variable v — & l n , the 'energy' by the non-relativistic formula 
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E = ^mv 2 , and the inverse temperature is (5 — mr 1 . This nonextensive gas 
has the special property that the temperature coincides with the mass of the 
'particles' considered. For nonzero a a similar formalism remains valid, just 
that the effective energies E become more complicated [22]. 

We can now evaluate all interesting thermodynamic properties of the system 
using the escort formalism. Regarding the invariant density (14) as an escort 
distribution P(&), the original distribution is given by 

= -Vl-$ 2 ~ P($) 1A? . (33) 

7T 



For the Tsallis entropy of the chaotic fields we obtain from eq. (24) and (33) 
SM = ~y - 2) = S q [P] (q = -l) (34) 

It is invariant under the transformation p — > P (a distinguished property of 
the entropic index q = ±1). 

For the generalized internal energy we obtain 
f 1 1 

U q [P] = j P(0)-m$ 2 = -m (q = -1) (35) 

-l 

and for the generalized free energy 

771 

F q = U q -TS q = T (3-vr 2 ) (q = -l). (36) 



All expectations formed with the invariant measure can be regarded as cor- 
responding to escort expectations within the more general type of thermo- 
dynamics that is relevant for our chaotically evolving scalar fields. Whereas 
ordinary matter is described by a statistical mechanics with q — 1, the chaotic 
fields generating dark energy are described by q — — 1. The role of dark energy 
is similar to that of a suitable thermodynamic potential in this more general 
type of statistical mechanics. In general, one can easily verify that the entropic 
indices q — +1 and q — — 1 are very distingushed cases: Only for these two 
cases the Tsallis entropy of the escort distribution is equal to the Tsallis en- 
tropy of the original distribution, for arbitrary distributions p(x). If our model 
of dark energy is correct, then both types of statistical mechanics are realized 
in the universe. 
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